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£NJ ■ Abstract 

We study the periodic Hartree-Fock model used for the description 
. of electrons in a crystal. The existence of a minimizer was previously 

shown by Catto, Le Bris and Lions {Ann. Inst. H. Poincare Anal. 
f"i . Non Lineaire 18 (2001), no. 6, 687-760). We prove in this paper that 

any minimizer is necessarily a projector and that it solves a certain 
nonlinear equation, similarly to the atomic case. In particular we show 
that the Fermi level is either empty or totally filled. 

1 Introduction 

On 

, The Hartree-Fock model is widely used to describe quantum electrons in 

usual matter. The purpose of this paper is to provide interesting properties 
CO ■ of the Hartree-Fock ground state in the periodic case, that is to say in a 

condensed matter setting. 
^ \ In Hartree-Fock theory [lOl [2] , the state of the electrons is represented 

by a so-called density matrix 7. This is a self-adjoint operator < 7 < 1 
acting on the physical space L 2 (M 3 ). When 7 has a finite trace, it models 
a finite number of electrons. An important example is given by a so-called 



Hartree-Fock state 

N 



l/ = ^2\<Pn)(V>n\ (!) 
n=l 

which is an orthogonal projector of rank and trace N, {<p n )n=i being an 
orthonormal basis of the range of 7. This example models N uncorrelated 
electrons, i.e. it corresponds to an A^-body wavefunction of the form Vl/ = 
if l A ■ ■ • A ifN in the antisymmetric product /\± L 2 (IR 3 ). The only correlation 
present in such a wavefunction is that of the Pauli principle, i.e. the fact 
that ^ is antisymmetric. 
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The Hartree-Fock energy is a nonlinear functional denned for any density 
matrix 7 and it has been largely studied in the literature [10\ [T3l [U O |3] . An 
important result due to Lieb [9] is that the minimization of this functional 
over all possible density matrices having trace N, gives the same minimum 
as when it is restricted to projectors of the form ([T]). Additionally, any 
minimizer is automatically a projector, 7 = Yln=i l^nK^nl) aim ^1) ■■■i l PN 
are the first eigenfunctions of a self-adjoint operator depending on 7 



a(H^) = {Ai < • • • < Ajv < • • • } U [0, 00). Notice that projectors of rank TV 
form the extremal points of the convex set of all density matrices of trace 
N. Although it may seem natural that a minimizer is always an extremal 
point of the variational convex set, this is indeed a non trivial result as the 
Hartree-Fock energy is itself not concave. All these properties are important 
for the device of efficient numerical methods [6j 0] . 

It was shown in [3] that necessarily Atv < A7V+1 which means that we 
can write the nonlinear equation satisfied by a minimizer in the form 



Throughout this article, we denote by \i the characteristic function of the 
set / C 1 and by Xi(A) the spectral projector on / of the self-adjoint 
operator A. 

This paper is devoted to the study of the periodic case, modeling an 
infinite quantum crystal. Our main goal is to prove that any Hartree-Fock 
minimizer satisfies properties similar to the molecular case described before. 
For the sake of simplicity, we assume that the crystal is simply the lattice 
1? and that there is only one nucleus of charge Z at each site of Z 3 . It is 
straightforward to generalize our result to any other periodic system. 

In the periodic case studied in the present paper, a state of the system is 
also described by a density matrix 7. But this is no more a trace-class oper- 
ator (it describes infinitely many electrons in the crystal) and it commutes 
with the translations of the lattice. As before, the orthogonal projectors, 
the extremal points of the set of all density matrices, will play a special role. 

In the periodic setting, there is also a periodic Hartree-Fock-like energy 
functional depending on 7, see Formula (jHJ) below. It was proved by Catto, 
Le Bris and Lions [7] that this energy admits a minimizer, but very few 
properties of this minimizer are known at present. We prove in this paper 
a result similar to [9j [3j. Namely we show that any minimizer 7 of the 
periodic Hartree-Fock energy is indeed always a projector and that it solves 
an equation of the form 



Hjlfn — \ n ip. 



(2) 



7 = X(-oc,x N ](Hy). 




(3) 
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with e 6 {0, 1}. The spectrum of is composed of bands and [i may be an 
eigenvalue (of infinite multiplicity, due to the invariance by translations of 
the lattice). Hence the nonlinear equation ([3]) solved by a minimizer cannot 
a priori be written in a simpler form like ([2]). 

In [Sj, a result similar to ([3]) was recently proved for a simpler model 
where the only non-convex term, called the exchange term, is neglected. In 
this case the proof is much simpler than in the general case treated in the 
present article. It is even possible to show that the spectrum of the self- 
adjoint operator is purely absolutely continuous, hence fj, cannot be an 
eigenvalue and one can take e = in ([3]). 

2 Model and main result 

For the sake of simplicity, we assume that the nuclei are placed on the 
lattice Z 3 and that they all have a charge +Z. We also discard the spin of 
the electrons. Our result can of course easily be generalized to any periodic 
system of spin- 1/2 particles. 

2.1 Notation 

Let fj be a Hilbert space. We denote by B(Sj) the space of bounded operators 
and by S(Sj) the set of (possibly unbounded) self-adjoint operators acting on 
$). In the whole paper, we denote by © p (fj) the Schatten class of operators 
Q acting on the Hilbert space Sj and having a finite p trace, i.e. such that 
tr^(|(3| p ) < oo. Note that ©i(fj) is the space of trace-class operators on 
Sj, and that ©2 (-ft) is the space of Hilbert-Schmidt operators on Sj. We 
note by L 2 (M) the space of square-integrable complex-valued functions on 
the Borel set M. Let us also introduce T = [— 1/2, 1/2) 3 the unit cell and 
r* = [—it, vr) 3 the first Brillouin zone of the lattice. 

In this paragraph, we introduce two functions G and W which we shall 
need throughtout the paper. They will respectively yield the so-called direct 
and exchange terms of periodic Hartree-Fock theory. We start by introduc- 
ing the Z 3 -periodic Green kernel of the Poisson interaction [TT], denoted by 
G and uniquely defined by 



and where the first equation holds in the distributional sense. The Fourier 
expansion of G is 




< 



\fcez 3 / 
G Z 3 -periodic, min G = 

R 3 



G(x) =h + 



fce27rZ 3 \{0} 
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with h = J G > 0. The electrostatic potential associated with a Z -periodic 
density p G Lj^M 3 ) n L 3 oc (IR 3 ) is the Z 3 -periodic function defined as 

(p * G)(x) := J G(x - y) p(y) dy. 

We also set for any Z 3 -periodic functions / and g 

D V> 9) ■= f J G(x - y) f(x) g(y)dx dy. 

Next we introduce the following function [7J 

W M = Enn-n> V,zeR 3 . (5) 
kez 3 |Z 1 

The function e lT, ' x W(r],x) is T-periodic with respect to x, when r] is fixed. 
So we can write W as a. Fourier series and obtain 

W(r,,x)=^e-^ £ (6) 

We have that, fixed x G T, W^(ry, x) — 4ir e ^ x is continuous in with respect 
to and that 

/ -iri-x\ 

lim f W{rj, x) - 47r ^p-J e" 7 '* + h = G ( x )- ( 7 ) 

2.2 The periodic Hartree-Fock functional 

The periodic Hartree-Fock functional was studied in [7J. 

The main object of interest will be the so-called (periodic) density matrix 
of the electrons. We define the translation operator acting on L^ oc (IR 3 ) 
as follows: Tfcu(x) = u(x — k) and introduce the following variational set of 
density matrices: 

V pei = | 7 € 5(L 2 (R 3 )) | < 7 < 1, Vfc G Z 3 , r fc7 = 7 r fc , 

tr L?(r) ((l - A^/\(l - A^ 2 ) <%<oo\. 

In the whole paper, we use the notation (^)^ e r* for the Bloch waves de- 
composition of a periodic operator A, see |15l 17]: 

A=-±fJ Azdt, Af€S(L}(T)), 
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L|(r) = ju € ^i 2 oc (^ 3 ) I T-fe-" = e~ lk <u, Vk G Z 3 } 

which corresponds to the decomposition in fibers L 2 (M 3 ) = f®„ d£L%(T). 

For any 7 G "P per (and almost every £ G F*), we denote by 7g(x,y) the 
integral kernel of 7^ G ©i(L|(T)). The density of 7 is then the non-negative 
Z 3 -periodic function of ^(M 3 ) n £ 3 oc (K 3 ) defined as 



Notice that for any 7 G V ] 



per 



/> 7 (x)dx = -—^3 ^ tr L 2 (r) (7 5 )^, 



i.e. this gives the number of electrons per unit cell. Later we shall add 
the constraint that the system is neutral and restrict to states 7 G V pei 
satisfying 

J p- 1 {x)dx = Z 

where Z is the charge of the only nucleus in each unit cell. 
The periodic Hartree-Fock functional is defined by 

£ (7) : = J tr L 2 (r) (~ A ?7 ^ ^3 - Z G{x)p 1 {x) dx 
r* 

+ lD(p 7 ,p 7 )-Ix(7,7) (8) 

for any 7 G V pe r- In the above formula, D(p 7 , /5 7 )/2 is called the direct term, 
while X(j, 7)/2 is called the exchange term. The latter is defined for any 
(3, 7 G Vpcr as 

*(/3,7) = ^e Jf^^'Jf d *dy P(Z,x,y)W(Z-e,x-y)W^yj (9) 



r*xr* rxr 



where VF is the function defined in ([5]). We remark that W{— rj, z) = 
W( V , z) = W(rj, -z), so X(J3, 7) = X(^J) = X( 7 , /3). 

In the whole paper, we use the convention that for a density matrix 

7 G P per 

1. \ 1 /" / /—j- /— t-\ d£ 



(2tt) 3 



which is well-defined by assumption. 
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In ([8]), we have considered that one pointwise nucleus is located at each 
site on Z 3 . We can also consider extended nuclei, in which case the last term 
of the first line of ([8]) is replaced by 

V per (x)p 1 (x) dx 

where 

V pei = Z T k m * G, 

m being a Cq°(R 3 ) nonnegative function such that L 3 m(x) = 1 and with 
a support small enough such that all the 77m have disjoint supports. For 
the sake of simplicity, we shall restrict ourselves to pointwise nuclei, the 
extension to the smeared nuclei case being straightforward. Of course the 
pointwise case corresponds to taking m = 5q. 




2.3 Existence and properties of minimizers 

The following was proved in [7], see Theorem 2.3 p. 698: 

Theorem 1 (Existence of minimizers [7]). Assume that Z > 0. The func- 
tional £ is well-defined and bounded from below onV vex . Additionally, there 
exists a minimizer 7 £ "P per of the minimization problem 

1= inf f( 7 ). (10) 

f T Pl = Z 

Notice in (|10p we could consider other constraints of the form J r p 7 = A 
but this does not make too much physical sense as the periodic system should 
be neutral in the thermodynamic limit. 



The purpose of this paper is to show that any minimizer 7 of (jlOD solves 
a specific nonlinear equation. Our main result is the following 

Theorem 2 (Self-consistent equation and the last shell). Assume that Z > 
and let 7 be a minimizer of (|10p . Then 7 solves the following nonlinear 
equation: 

7 = X{-oo,n){H^) + ex{^}(^ 7 ), 

{H 7 ) ( = - A e - ZG + Pl * G - (2tt)- 3 j ^ W{t-' - e, x - y)l^(x, y) d£', 

(11) 

where e G {0, 1} and /i € M is a Lagrange multiplier due to the charge 
constraint J r P-y = Z. 

Remark 1. In the periodic case, for the minimizer 7 to be a projector, Z 
need not be an integer as in the atomic setting. 
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A result similar to Theorem [2] was proved for the reduced model in [5, 
Theorem 1]. The reduced model consists in neglecting the exchange term, 
i.e. the last term of ([8]). In this case, the Hartree-Fock energy becomes 
convex and it admits a unique minimizer. The proof that it satisfies an 
equation similar to (lllj) is then much easier because the only term which 
depends on 7 in the operator is p~ * G when the exchange energy is 
neglected. Using that p 7 * G is just a periodic function (hence its Bloch 
decomposition does not depend on £), one obtains from a result of Thomas 
[T7] that the spectrum of is purely absolutely continuous. In particular 
there cannot be any eigenvalue, hence we can take e = in (llip . 

When the exchange term is kept in the model as in the present paper, 
the situation is much more complicated. The main difficulty is that the 
Bloch decomposition of the last term of depends in a non trivial way of 
£. In particular, we do not know if the spectrum of is purely absolutely 
continuous. It is in principle possible that \x is an eigenvalue (of infinite 
multiplicity) of H~ in (jlip . However, we are able to prove that any mini- 
mizer 7 is automatically a projector. Theorem [2] even states that either the 
minimizer 7 does not contain the eigenspace corresponding to the eigenvalue 
fj, (e = 0) or it fills it completely (e = 1). The fact that 7 is a projector is 
an important property which can be used as a basis for the construction of 
models for crystals with local defects, as this was done in [5]. 

The rest of the paper is devoted to the proof of Theorem [2j 



Lemma 1 (Properties of W). We recall that W is defined in ([5]). We have 
for all (£, x) E V* x V 

x) = 4ir^- + e-^{G{x) - h) + e^/fr x) (12) 



3 Proof of Theorem [2] 



Step 1. Properties of W . 



We start with the following useful 



where f satisfies 



V£,f E r* 



L°°(r) 



<cie-e 



(13) 



Vxer, /(o,x) = o. 

Similarly, we have for all (£,£', 2, y) E (r*) 2 x T 2 



(14) 





<i 




+ g-^-e'M* 



v)G{x-y) + g{i-i\x-y) (15) 
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where g G L°° ((r*) 2 x T 2 ) safe/ies /or aM G T* 

sup |<?(6 - -5(6 <c|6-£i + (16) 

x,yeT 2 

In (|15p . we have used the notation \x\oo = sup{|xj|, i = 1,2,3}. 

Proof. The proof is contained in [7\ p. 744] and we only sketch it for the 
convenience of the reader. A simple calculation gives 

mGZ3\{0} IS III 

Hence we obtain for e small enough 

||/(e,-)-/(e',-)|li 2((1+e)r) <^-e'i 2 . 

Next one uses that x \— ► /(£, cc) — /(£', x) is harmonic to infer 

|/& •) - /(£', oiLocfr) < c •) - /^', oILi^d ^ ^ - e'l- 

The proof is the same for g. □ 
In view of the previous result, it is natural to introduce 

X G ( 7 , 7 ):=// dtdZ'ff dx dy G(x - y)e~^-^< x -^^(x, y)^ (x, y) 

dxdyG(x -y)\^(x,y)\ 2 (17) 

rxr 

where we have defined j^(x, y) = e~ l ^' x ^^(x, y)e l ^' y . Remark that p 7 = P j. 
Definition (|17p is quite natural as it is an exchange term taking the same 
form as in usual Hartree-Fock theory [10] . For all 7 G "Pper we have 

h(x,y)\ 2 < Pl {x) Pl {y), (18) 

which is a consequence of 7^ > 0, see, e.g., [7, p. 746]. The same inequality 
holds with 7 replaced by 7. We infer that 

*g(7,7) <D{ Pl , Pl ). (19) 

Indeed 7 1— ► -D(p 7 , p 7 ) — Xq^, 7) is easily seen to be continuous for the weak 
topology by Fatou's Lemma, and this can be used to prove the existence of 
a minimizer. 
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Step 2. Regularity of the spectral decomposition of 



In the rest of the proof, we consider a minimizer 7 for (jlOD . Such a minimizer 
is known to exist by Theorem [H proved in [7]. We know that it has a finite 
kinetic energy since it belongs to Vp er : 



JL 



< 00. 



We recall the useful inequality Eq. (4.42)] 



tr 



r* 



(27r) J 



> 



(20) 



^6iri (r)cLi (r)n^ er (r) 



IS 



which proves that the map 7 £ "P pcr 
continuous^- 

We investigate the regularity of the eigenvalues and eigenvectors of H~ r 
For simplicity, we study on the fixed Hilbert space Lp er (r). This means 
we introduce the unitary operator defined in each Bloch sector by 

2 ^ > L|(r) 

u(x) 1 — ^ e l *' a; «(a;) 



We shall use the convention that when T£ is an operator on L|(r), then 
:= UfT^U^. Let us denote by X the exchange term defined by its kernel 

X^x, y) = (2vr)- 3 J ^ W(e -£,x- y)^(x, y) df. 



Then we get using ([T5 



4. ^ 



<i 



X^(x,y) = U^X^(x,y) = 

+ G(x - y) + <?(£ - x - y) 
where we have introduced g(rj,x) = e iri ' x g(rj,x). 



^liTtm-(x-y) 

|^-^-2vrm| 2 



(2vr) s 



(21) 



Lemma 2. The family (X^)^ e r is bounded in B(L^ er (T)) and it is Holder: 



V0 < p < 1, 



x 6 - X 6 



B(i?«(T)) 



<c P i£i-6r. 



1 This property can be used to properly define the term D{p 1 ,p 1 ) and see that it is 
continuous for the topology of /C. 
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Proof. We have for all u,v G L 2 er (T) 



<4vr £ 

<1 



|£' -£i - 2vrm| 2 |£' -£ 2 - 2vrm| 2 



X \{^'U2 7 rmU,U 27 TmV)\ 7^33 



+ 11 dxdy \g(£' - £ l} x - y) - <?(£' - fr,x - y)\ x 
J Jrxr Jr* 



x \%'(x,y)\ \u{y)\ \v(x)\ 



(2tt) e 



We recall that < 7 < 1 which means that < 7^ < 1 for all £ G T*. Using 
(fl~6l) we obtain 



<C|6-6I 



rxr 



+ 4 ' £ /r. 

<1 1 

As < 7; < 1 for all £ G T, we have 



(2vr) 3 



-2vrm| 2 |£ / -£ 2 -27rm| 2 



(2tt) e 



(2vr) 3 dxdy\^,(x,y)\ 2 = / tr £ 2 cr(r) [(% 

t/r* 1/ 1/ r x r i/r* 



(2vr) s 



This easily proves the Holder regularity. 

The proof that (X^)^ £ r is bounded in £>(L 2 er (T)) is essentially the same 
for all the terms which we have treated above. We just estimate the term 
involving G(x — y). We notice that 

2 



G(x - y)%'(x,y) 



rxr 



(2vr) 3 
dx dy G(x — y) 2 



e2(L| er (D) 



G(x -y) 2 \^(x,y)\ z dxdy < II G{x - yf p^ f (x)p J (y)dx dy 



rxr 



rxr 



by (fl8j) . The last term is well defined (and it is independent of £). To see 
that, we notice that for all x, y G T 

1 



G(x - y) 



E 



meZ 3 , I <1 



\x — y + m\ 



h(x - y) 
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where h(x — y) is bounded on T x T. We treat for instance the term with 
m = 0, the argument being the same for the others. Fix some non negative 
cut-off function x which equals 1 on T and vanishes outside (1 + e)T for e 
small enough. We have by Hardy's inequality \x\~ 2 < 4(— A) and using the 
periodicity of p 7 together with the formula J r p 7 = Z 



IL 



P^)pM dxdy < [ dx f dy X(x) 2 Pl (x)p,(y) 



rxr \x-y\ 2 Jm? J r \x-y\ 2 

<AZ dx\V(x^)(x)\ 2 <C [Z + / dx\V^(x)f 

Arguing similarly for the other terms and using (|20l) we get 



G(x - y)^'{x,y) 



(2vr) 2 



e 2 (^ er (r)) 



Kcfz + Jd^tr^ [(-A)^ 



which is independent of £. □ 
Recall that the operator (Hy)^ = U^(Hy)^U^ acting on L 2 er (T) reads 

(# 7 ) e = -^-iS.V + \t£-ZG + p 7 *G-Xt (22) 

where of course —A and V are respectively the periodic Laplacian and the 
periodic gradient acting on Lp er (r). We denote by A&(£) the eigenvalues of 
(Hy)^ which are the same as that of (H~)£ since is unitary. We may 
assume that the Afe(£) are in nondecreasing order: Ai(£) < A2 (£)<••• . 

Lemma 3. The eigenvalues Afc(£) tend to +00 when k — ► +00, uniformly 
in £ G r*. 

Proof. We have proved in Lemma[2]that is a bounded family of operators 
in B(L 2 er (T)). Also we have p 7 > and G > 0. Hence the following holds 
on L 2 pei (T) 

(H 7 )t > -A/2 - i£ • V + |£| 2 /2 -ZG-C. 
As G £ £ 2 (T), it suffices to apply [13 Lemma A-2]. □ 

Lemma 4. Consider f2 an open set of T* and let K be a compact set in C 
such that mi^ € Qd(K,a(Hy)^) > e > 0. Then we have that 

1. (1 — A) ^(if 7 )^ — z^j is bounded on L 2 CI (T), uniformly with respect 
to £ £ Q and z G if; 
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2. the map £ i-> (1 - A) ^(H^ - z 
respect to £ G £1, uniformly in z € K 



G ^(Lp er (r)) is Holder with 



Proof. By Lemma [31 we know that there exists a real number zq > such 
that d( — zq, a(Hy)A > e for all £ G T* and some e > 0. We choose a 
constant c> zq large enough. We have 



(tf 7 ) € + c) (c- A)" 1 = l/2 + i? c 



where 



V 



Rc ■= x + |£| 2 /2(c - A)- 1 + { Pl *G-ZG- X{) (c - A)" 1 . 

c — A 



We have 



V 



c - A 



6(-£| er (r)) 



+ |||£| 2 /2(c-A)- 1 | 



< 



B(^ er (r)) - j + 



and, as is uniformly bounded in £>(Lp er (T)) by Lemma [2] 



C 



IWc-A) || B(L , cr(r)) < — 
Next we can use that G G Lp er (r) to get that 



|G(x)(c- A)" 1 ! 



S2(ig er (r)) 



^(dIK^h 2 )- 1 "-' 



< 



c 

„l/4 ' 



Similarly we have G G L per (r) and p 7 G Lp er (T), hence p 7 * G G Lp er (r) 
and the same estimate furnishes the bound 



\\( Pl *G)(x)(c-A) le.CL^cr)) 
Eventually we have proved that 



< 



C 

-1/4 ' 



l»(i?«(r)) 



< 



C 



-.1/4 ' 



This shows that for c large enough, the operator ((-ff 7 )g + c)(c — A) 1 is 
invertible and that its inverse is bounded: 



(23) 



(c-A) (fT 7 ) € + c GS(^ er (r)). 



The last step is to use the formula 



-i 



-i 



JET 7 ) e + c) + (z + c) I (H 7 ) 5 + c 



-i 



-i 
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which, when inserted in (|23[) yields that (1 — A)((i/ 7 )g — z)~ l is bounded 
on Lp er (r), uniformly with respect to £ G and z £ K. 

For the second point of the Lemma, we use the resolvent formula 

(1-A) 



1 - A) - - (# 7 ) e ) ((tf 7 ) e 



2 



-i 



(i-A) (# 7 ) c - * - O • v + \e - \C\ 2 - x ( + X; 



x(l-A)- 1 (l-A)((i? 7 V 



-i 



Next we can use that £ t—* is Holder for the norm of B(Lp er (T)) as shown 
in LemmaO and that V(l — A) -1 is a bounded operator on Lp Cr (T). □ 

Lemma 5. The eigenvalues Afc(£) are Holder with respect to £ G T. 

Proof. As before we take c large enough such that d( — c,a(Hy)A > e 
for all £ G T* and some e > 0. For every £ G r*, the spectrum of the 
self-adjoint operator := ((H^)^ + c) _1 is composed of the eigenvalues 
MfcCC) := (Afc(£) + c) -1 . By Lemma [U £ i— > i?(£) is a Holder family in 
0(Z,£ er (r)). By the usual min-max Courant-Fisher formula, we have 

iM*(e)-/i*(Oi < \\r(o - 

hence for any k > 1, £ € T* i— ► /ifc(£) is a Holder function. This a fortiori 
proves the same property for the eigenvalues Afc(-). □ 

Lemma 6. Let f2 6e an open subset of T* and I = (a, b) an interval of R 
smc/i that cr(iJ 7 )^ n {a, 6} = /or a// T/ien i/ie map 

is Holder. In particular, we can find an orthonormal basis («i(£), ...,uk{0) 
of the range of xi(Hy)z such that £ G SI ^ u&(£) G i^p Cr (r) is Holder with 
respect to £ G $7. 

Proof. This is a simple application of Cauchy's formula [8] 

Xi(H 7 )t = ^-j c dz - z) 1 (24) 

where C is a smooth curve in C enclosing the interval I and intersecting the 
real axis at a and b only. □ 
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Step 3. Variational property of a minimizer 7. 

We follow a well-known argument pQ, [2] (see [121 Section 4] for a very similar 
setting) and consider a fixed state 7' € Vp eT := {7 6 "Pper | J r p 7 = Z}. As 
P^gj. is convex, we have (1 — t)j + £7' S "P^f er for all t G [0, 1]. Hence, since 
7 is a minimizer 

v te[0 ,i], Eh±M^m> . , 25) 

Expanding, we obtain 

b W (W-7){)^j + 

+ ^|z?(py„ 7 ,py_ 7 )-X(7 , -7,7 , -7)| >0 (26) 

where -ff 7 is the mean-field operator defined in (jlip . Taking i = 0, we obtain 
that 7 is also a minimizer of the linearized functional at 7, 

(2vr) 3 ' 



r* 



i ^ t r L|(r) ((-Ky)^) 



r* 

on the convex set Vp er . For any £ € T* we can choose {<£&(£, ■)}& a basis of 
eigenfunctions of (H^)^ with eigenvalues Afc(£), such that 

(H^ = J2^k(0\fk(0)M0\- (27) 
fe>i 

We know from Lemma [5] that each Afc(£) is Holder with respect to £. Let us 
introduce like in Appendix] the function 

*!>i 

The function C is nondecreasing on R. The operator i/ 7 being bounded 
from below, we have C = on (— 00, inf Afe(T*)). Also we know that 
limre-^oo C(k) = +00 by Lemma [3l We deduce that there exists a \i S R 
and a periodic operator (5 € Pper such that 

lim C(k) <Z< lim C(k) 

k— >/i~ re— 

and 

7 = X(-oo )At )(^ 7 ) + ^ 
where < 5 < 1 and Ran(5) C ker(ff 7 — /i). The proof of this fact was 
given in the Appendix of [5j. When [i is not an eigenvalue of -ff 7 , i.e. 
|{£ G T* I 3k, Afc(£) = fi}\ = 0, the proof of Theorem [2] is finished, as nec- 
essarily 5 = 0. What rests to prove is that if [i is an eigenvalue of -ff 7 , then 
S = or 5 = X{f,}(H~t)- 
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Step 4. The Fermi level is either empty or totally filled. 

Now we argue by contradiction and assume that \x is an eigenvalue of H^: 

liter* \3k>i, A fc (0 = ^}|^o. 

We also assume that 5^0 and 5 ^ X{ix}(H~t)- The following lemma will be 
a key result to construct perturbations of 5. 

Lemma 7. Assume that \i is an eigenvalue of H^, that <5 7^ and 5 ^ 
X{fi}(H~/)- Then there exists a constant e > 0, a Borel set ACT* with 
\A\ ^ and two continuous functions ^6ii-> u(£) E Hp eT (T) and £ E j4 i— > 
u'(£) G i?p er (r) suc/t i/wrf 

V£€A, u(0V(0 Gker((# 7 ) c - M ), 

l«(0ll^ ef( r) = lh'(0|L2 ef(r) = l 
and, denoting ip(£) = U?u(£) and <p'(() = U^u'(^), 

0<5s + t\<p(Z))(<p(Z)\-t'\<p'(0)(<p'(0\ <1 (28) 
on L|(r), /or all £ £ A and all t,t' £ [0, e). 

Proof. Assuming the eigenvalues are in nondecreasing order and using that 
Afe(£) ~~ * 00 as ^ — > oo, uniformly in £ and Lemma [5l we deduce that there 
exists fei, m, e' > and a subset ^4 of T* with |A| 7^ such that 

A fcl -i(0 < A* - e' < A fel (0 = m = A fel+m _i(0 < /i + e' < A fe „ +m (0, (29) 

< e' < tr L 2 (r) (n 575 n ? ) < m - e' (30) 

for all £ G ^4 and where we have introduced 11^, the orthogonal projector on 
ker((i? 7 ) € -/i) in Lf (r). As < 75 < 1 for all £ G T* and n^n^ = c%, ([30]) 
means exactly that 5^ is neither nor 1 when £ £ A. 

Now we choose an adequate basis of the range of LT^ . By Lemma [SJ we 
know that the eigenvalues A&(£) of (H^)^ depend continuously on £. Hence 
(f2"9"|) is satisfied on the set A. As \A\ 7^ 0, we deduce that the interior of A is 
not empty. Hence, decreasing A if necessary, we can assume without any loss 
of generality that A C where f2 is an open set on which (|29|) holds tru^l. 
By Lemma El we can choose u^i^), nfc 1+m _i(£) an orthonormal basis in 
Lp er (r) of ker((H^)^ — fi) for all £ in ft such that £ 1— ► n/%(£) is a Holder 
map in i7p Cr (r) on $7. Let us recall that (i? 7 )^ := U^(Hj)^U^ where C/g : 
Lp er (r) — > L|(r) is the unitary operator which acts as a multiplication by 

2 However we cannot assume a similar property for (|30p because we have no information 
on the regularity of 7^ on the range of Ilf . 
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the function e l ^' x . Of course the functions <pk(£) '■= U^Uk(£,) are eigenvectors 
of (-ff 7 )^ with eigenvalue //, by definition. 

Next, for any £ £ A, we may introduce the m x m matrix of n^II^ 
in the basis {MO^Jk?^ 

{M^ij = (^i-i+iCO'^Vfei-i+iCO)^^)- 

As £ i— > £ 5(C m ) is a bounded measurable function, Lusin's theorem 
|16| tells us that it is continuous on a Borel set B C f2 with |f2 \ £?| as 
small as we want. Hence, replacing A by An B if necessary, we may assume 
that £ I— ► Mg G 5(C m ) is continuous on A. The eigenvalues < Wi(£) < 
• • • < n m {^) < 1 of the matrix Me are then continuous functions on A. Our 
assumption is that on A, Me is not the zero matrix and is not the identity 
matrix. 

First we treat the case when Me has an eigenvalue which is not equal to 
and not equal to 1. This means decreasing A one more time if necessary, 
we may assume by continuity that for all £ G A 

rik(Q < e < n k+1 (0 < ■ < rik+JxiQ < 1 - e < n k+Jl+1 (£) 

with the convention that no = — oo and n m +i = oo. By Cauchy's formula 
(pill the projector on ®^=i ker(Mg — n^+j) depends continuously on £ G 
A, see [8]. Hence we can find a continuous function £ G ^4 i— > v(£) G 
0/=i ker(M 5 - n fc+i ) C C m such that 

Vt€(-e,e), < Mf + < 1- (31) 

and |v(£)| = 1 for all £ G A We may then introduce 

m 

^) = E"(OiW 1 -l+i(Oei£(T) and U (0=%(0€^(r). 
i=i 

Notice that the function £ G ^4 ^ u(£) G iJp er (r) is continuous by the 
definition of v(£) and To get ([28]) . we just use <p' := 93 and u' := u. 

The other case is when Mc is always a projector for almost every £ G 
A By continuity of £ 1— > tr(M^) we may, decreasing ^4 one more time if 
necessary, assume that it has a constant rank < J2 < m. This means we 
have on A 

= ni (£) = ••• = "m-J'(0 < n m-j 2 +i(£) = • • • = n m (0 = 1. 

As before by Cauchy's formula (|24|) the projectors on ker(Mg) and on ker(M^— 
1) depend continuously on £ G A. Hence we can find two continuous func- 
tions £ G A i-> G ker(M f ) and £ G A h-» G ker(M f - 1) such 
that 

Vt, t' G [0, e), < M e + MO* " ^'(6^(0* < 1 
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and |t> (£)| = = 1 f° r an £ £ A. It now rests to define 

m 

3=1 
m 

Equation (J2SJ holds with = ^(C) and «'(£) = f7$ ¥>'(£)• □ 

Lemma 8. Assume that p is an eigenvalue of H^, that <5 7^ and 5 ^ 
X{ At }(^ 7 ), and let be A and <£>(£), <p'(£) as in Lemma\7\ Then for all rj,rf G 
L°°(A,]R + ) such that f A 'H = f A rj', we have 

D(p R ,p R )-X(R,R)>0 (32) 

where R is the periodic operator defined by 

= v(0M0)(v(0\ - v'(0\<p'(0W(0\- 

Proof. Let 77,7/ G L°°( J 4,R + ) be such that J A (j] — rf) = and R as defined 
above. First we remark that for < (l^ll^oo + < e -> we have 7 + 1R G 

7^f er by construction of </?(£) and ¥>'(£)• Next we use Equation (|26j) with 
7' = 7 + ii?. We compute 

d£ /• ..,.,„ . ... .... d£ 



J tr L 2 (r) [(fl" 7 ) { ^ = J r,(Z){(H 7 )z<p(Z)M0) 
r* a 



L «( r ) (2vr) s 



l/(0<(fl7W(0^(0>^CT)^8 



(2tt) s 

A 

hence the first order term of (|26|) vanishes and the result follows. □ 

The last step is to construct functions 77, rj G L°°(j4, M + ) such that J A (rj— 
7/) = and D(p R ,p R ) — X(R,R) < with R defined as above. This will 
contradict (|32p and finish the proof. Indeed, we even prove that 

inf (D(p R , p R ) - X(R, R)) = -00. 

jj,77'GL 00 ( J 4,R+), 

f A (v-v')=o 

The idea of the proof is somewhat similar to that of the atomic case [9, 
[H El [21 [12]. The role of the perturbation R is to transfer some charge 
from the eigenvector </>(£) to the eigenvector y>'(£) within the last level p. 
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The originality of the periodic case studied in the present paper is that the 
transfer needs to be done between two different Bloch sectors L? (T) and 
L| 2 (r). This means essentially we want to take 77 = 5^ and rf = 5^ 2 in 
Lemma El 

As \A\ 7^ 0, we may find two points £1 and £2 in A such that \A n 
Bfa, A) I ^ and |A n A)| / for all A > 0. Here £(£, A) denotes the 
ball of radius A centered at £. We may also assume that < |£i — £2! < 1/4 
in such a way that £ — £' G T* if £ G -B(£i, A), £' G -B(£2, A) and A is small 
enough. Next we define 



1 



Vx 



A-lB(a.A) 

|A n S(a,A)|' 



1 



Anff(&,A) 

|Ans(£ 2 ,A)| 



in such a way that t?a — 1 % and rj' x — weakly as A — > 0. We denote 
by i? A the associated periodic operator, the Bloch decomposition of which 
is given by 

We also introduce as usual the family of operators acting on the fixed space 

-kper(T) 

Next using (fl2j) . we get 



(27rf(Z)(^,p^)-X(i? A , J R A - 



/ (A)+/ 1 (A)+/ 2 (A) 



where 



/o(A) = (2kY{D{ Prx , Prx )-X g {R\R x - 



+ h 



rxr 



dx dy 
dx dy 



77 A (£)«(£, *)«(£, y) - r/ A (£)u' (£, x)u'(£ ; j/) J d£ 



rxr 



ixi 



x hA(£)w(£,x)u(£,y) -r? A (£)u'(£,;c)u'(£,y) x 



+ 8vr 



d£d£ 



x {Vx(t>^,x)u(e,y)- V ' x (e)u>(t>,x)u'(e,y) 

,vx(0v' x (e^ 2 



AxA 



-4tt 



ie-ci 2 

, Vx(£)vx(£' 



AxA 



ie-£ 



<«(0,« / (O>xg« t r) 

M6XO>L| er( r) 
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/ 2 (A) = -4tt // ^^ V \f W ?P |<^(0y(O>^.cr) 

To finish the proof, it rests to show the following 
Lemma 9. We have 

lim 7 (A) = J € R, lim IUX) = lim J 2 (A) = -oo. 

A^O A^O A^O 

Proof. Using the important property that the maps £ £ T* ^ u(£) € 
i7p er (r) and £ G T* G i?p er (r) are continuous by construction, 

we see that 

lim/ (A) = (^f J/ dxdyG{x - y) (Z(x,x)Z(y,y) - \Z(x,y)\ 2 ) 
rxr 

+ m II dXdy (j^ 1 ~^ x ~ v ^ + -^l 2 ) u (^ 1 ' a; ) u (^ 1 'f) x 
rxr 

xu'(£ 2 ,x)u'(£ 2 ,y) + /i ^ cfedy|Z(x,y)| 2 



rxr 

with 



Z(x,y) := u(£i,x)u(£i,y) - u'(£ 2 , x)u'(£ 2 , y). 

We have also used that f(0,x) = 0. Finally, we use again the continuity of 
£ i ^ u(£) in Lp Cr (r) to infer that 

<«(O,«(0>Lg er (T) = l+OA- + 0(l) 

when G B(£\, A). Hence, for A small enough 

M A)^//«'?M S -! 

AxA 

which proves that limA^o A (A) = — oo. The argument is the same for I 2 (A). 

□ 

Acknowledgement. We would like to thank Eric Sere and Isabelle Catto 
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